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Abstract 

We propose a model for describing the growth on 
an untreated tumor, which is characterized in a sim- 
ple way by a minimal number of parameters with 
a well-defined physical interpretation. The model is 
motivated by invoking the Master Equation and the 
Principle of Detailed Balance in the present context, 
and it is easily generalizable to include the effects of 
various types of therapies. In the simplest version 
that we consider here, it leads to a linear equation 
that describes the population growth in a dynamic 
environment, for which a complete solution can be 
given in terms of the integral of the growth rate. The 
essential features of the general solution for this case 
are illustrated with a few examples. 
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1 Introduction 

In the last few decades, the mathematical mod- 
eling of tumor growth as a function of time has 
been based mainly in applications of the Gompertz 
equation|T]ElEl and a modified version of it^, power 
law equations^ and generahzations of the lo- 
gistic equation|Bl [?]■ These mathematical models 
helped to provide an understanding of tumor growth 
as a more organized phenomenon than previously 
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thought. In particular, solutions of the Gompertz 
and the modified logistic equations exhibit an S- 
shaped form which represents inhibition of growth to 
an asymptotic limit. In addition, their applications 
to experimental data, with the help of additional dif- 
ferential equations to model therapy, has served as a 
guide to improve the effectiveness of treatment. 

However, a preference for one model or the other 
has been based exclusively on their adequacy to fit 
experimental data. Similarly, some of the models 
consist on generalizing previous ones by modifying a 
parameter in a way that is not motivated by any fun- 
damental principle, but again to fit some data. For 
example, the so-called Generalized Logistic model is 
the result of merely modifying the logistic equation 
by inserting an arbitrary power which, in order to fit 
the data, becomes a noninteger number. 

While this may be appropriate for particular pur- 
poses, and for applications to systems that are not 
governed by the laws of Nature, the fact that the 
models are neither motivated nor based on funda- 
mental physical principles prevents us from ascribing 
a physical meaning to the parameters that appear 
in such models. This makes it more difficult, if not 
impossible, to establish a connection between the ob- 
servable effects that can be described on the basis 
of phcnomcnological models, with more fundamental 
explanations and understandings of the mechanism 
of growth in biological systems, which must be based 
on detailed microscopic dynamics. 

The model is motivated by invoking the Master 
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Equation and the Principle of Detailed Balance in the 
present context. In the simplest version that we con- 
sider here, it leads to a linear equation that describes 
the population growth in a dynamic environment, 

The main objective of this article is to present a 
model that is inspired by fundamental physical prin- 
ciples, such as the Master Equation and the Princi- 
ple of Detailed Balance, and is partly motivated by 
analogous equations that apply to a variety of phys- 
ical systems. The model is based on a linear dif- 
ferential equation with a time dependent coupling. 
It shares the same asymptotic behavior of the GL 
model, giving an S-shaped function for tumor growth 
that vanishes at t = and eventually reaches an 
asymptotic value at large t. However, the model pro- 
posed here contains a minimal number of parame- 
ters, which have a concrete and well-defined meaning 
and are in principle determined and calculable if the 
interactions that govern the underlying microscopic 
mechanism of growth are known. By the same token, 
the use of this model should in turn shed light on 
such mechanisms, thereby providing a firm footing 
for pursuing and extending such approaches. 

The rest of this paper is organized as follows. In 
Sectional for the purpose of comparing our proposed 
differential equation with the previously proposed 
models, we briefly recall the basic aspects and short- 
comings of them. In Section Owe discuss the theo- 
retical framework that motivates the model that we 
propose, the assumptions and idealizations involved, 
and the interpretation of the parameters that appear 
in it. In Section^ the general solution to the equa- 
tion is displayed, and it is illustrated by considering 
various specific examples that can be of practical use. 
Finally, Section [3 contains our conclusions. 

2 Previous models 

We will denote by f{t) the growth function, and for 
simplicity of the notation omit the argument t except 
when necessary to avoid confusion. We envisage / to 
be the number of cells at a particular instant of time. 
The Gompertz equation can be written in the form 



4f 
dt 
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where a and P are two parameters. The solution is 
given by 

fit) = fie^ , (2) 



with 
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where fo is the initial value of /. The asymptotic 
value is determined as 



(4) 



A particular feature of Eq. ^ is that it is not de- 
fined for / = 0, and therefore it does not allow the 
initial value fo = 0, as Eq. (|3J) reveals, thus ruling 
out its application to data with rather small initial 
population value. 

A way to overcome this difficulty is to modify Eq. 
([T]) by introducing an arbitrary parameter e > such 
that the new differential equation reads 0] 



dt 



af 



In 
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which has a non trivial solution given by 



fit) = /3e 



where 



G'= — 



1 

eat 
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While this procedure overcomes the above-mentioned 
difficulty with the Gompertz solution, a physical mo- 
tivation or justification is lacking. 

Another approach|Sj has been to consider the 
power law differential equation 



dt 



For /3 ^ 1 it has the solution 

fit) = [/3i(l ~a) + /J 



l-Ql 



(8) 



(9) 



which leads to linear or exponential growth for a = 
or a = 1, respectively. Certainly, power law growth is 
unconstrained, and its behavior is radically different 
from that given by Eq. 
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Another way to obtain asymptotic behavior in tu- 
mor growth is to modify the logistic equation by in- 
troducing an arbitrary power e > 0, leading to the 
differential equation 



dt 



1 - 



whose solution is given by 



where 
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fo 
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fo 
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Eq. (|ll|l reproduces the Gompcrtz function for e = 0, 
but it also has the limitation that the differential 
equation Eq. H1U|) does not allow a zero value at f = 0. 
Therefore, the solution of Eq. H1U|I can be consid- 
ered to be the most general function that possesses 
asymptotic behavior and a shape consistent with ex- 
perimental data. However, the fact that Eq. (|10|l 
can be neither motivated nor understood on the ba- 
sis of some of fundamental physical principles is in 
our opinion an important limitation to further un- 
derstanding the mechanisms that drive the growth in 
these biological systems, on the basis of such equa- 
tions. The main reason for this lack of insight is due 
to the fact that Eq. H1Q(I is the result of merely mod- 
ifying the logistic equation by inserting an arbitrary 
power, that in order to fit the data becomes a nonin- 
teger number. 

Therefore, without any other justification it is not 
possible in this and the previously discussed models 
to ascribe a deeper physical significance to the pa- 
rameters of the model. 

3 The Model 

Our starting point is the equation that describes the 
growth of a population that is sustained by an envi- 
ronment. We assume that in such situations the pop- 
ulation grows up to a certain saturation limit fs , and 
that the environment is large enough such that it is 



not affected by the population itself. Under such con- 
ditions, we assert that the rate of change of the pop- 
ulation is proportional to the difference between ac- 
tual value of the population and its saturation limit. 
Therefore, 



/ = -7(/ - fs) 



(13) 



This equation is reminiscent of Newton's cooling law 
which states that the rate of change of the temper- 
ature of a system is proportional to the deviation of 
the system's temperature from the temperature of its 
environment. In our context it is possible and use- 
ful to give a motivation and justification in terms of 
more basic principles as follows. 

3.1 The Master Equation 

The problem of the time evolution of the population 
of a given specie appears in many physical contexts. 
For example, in the astrophysical context of the Early 
Universe, one analogous problem is the determina- 
tion of the abundance of the various atomic elements 
and how they form'S'. The nucleosynthesis processes 
in stars are examples of similar phenomena which, 
among other things, explain the generation of energy 
in the Sun. 

In one way or another, the two basic principles that 
guide the development of a population are the Master 
Equation and the Principle of Detailed BalancepillO|. 
The master equation takes the form 



(14) 



where W depends on / itself and the other variables 
that describe the rest of the system with which the 
population can interact. W is decomposed into a 
series of terms, each of which represents the contri- 
bution due to a particular process that causes the 
population to change. The principle of detailed bal- 
ance states that there is a precise relation between 
the so-called direct process and its inverse. 

For example, let us consider a process in which only 
one cell participates and let us denote such process 
in symbols by 

(15) 
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where (j) stands for a member of the population (a 
cell) and X stands for a different object. In the di- 
rect process, indicated by the right-pointing arrow, 
a cell (j) disappears into X, while in the inverse pro- 
cess, indicated by the left-pointing arrow, the reverse 
is true. Then, denoting by jd and ji the rates for 
the direct and inverse processes, respectively, their 
contribution to W is written in the form 

Wi = -7(/ - /.) , (16) 

where 

7 = 7<i - 7i • (17) 

Similar equations also describe the kinetic ap- 
proach to equilibrium of systems that are put in con- 
tact with a reservoir. In such cases, which are gov- 
erned by physical kinetic equations, the principle of 
detailed balance implies a fundamental relation 



where AE and T are identified with the change in 
energy of the system and the temperature of the en- 
vironment, respectively. 

In our case, in principle both "fd,i could be calcu- 
lated if the interaction between the cells with their 
surroundings and among themselves were known. 
However, we have at present no such theory of these 
interactions. Thus, we leave 7 as an unknown param- 
eter with the property that is a positive quantity. 

The procedure outlined above for the case of single 
cell processes can be generalized to more complicated 
ones. For example, consider the processes in which 
two cells participate, which we denote in symbols by 

(j)(t)^X. (19) 

Because the direct process involves two cells, its rate 
is proportional to By the same reasoning that 
lead us to write Eq. , the contribution from these 
processes to W is of the form 

W2 = -7'(/' - /.') , (20) 

where 7' characterizes the rate for the process to oc- 
cur. 

As a typical rule in those contexts in which these 
equations have already been applied, the processes in 



which more than two members participate are rare 
and not important. Therefore, we are tempted to 
state that the master equation 

| = -7(/-/,s)-7'(/^-/.^), (21) 

is a good starting point for further exploration of 
these ideas in the present context as well. 

In the present paper, we will restrict ourselves to 
the linear term only, as written in Eq. H13|) . The as- 
sumption behind this approximation is that the pro- 
cess in which the cells participate in pairs are rare 
compared to those in which only one cell participates. 
Should this linear approximation prove to be inad- 
equate, it could indicate that the pair interactions 
are important and the quadratic terms in Eq. (|?T|l 
should be taken into account. Overall, this approach 
provides a framework for carrying a systematic anal- 
ysis, based on incremental approximations, on a firm 
footing and in an organized fashion. 

4 Solution 

4.1 Static environment 

When 7 is a constant, Eq. has the simple solution 

fit) - /. [1 - e-^*] + foe--" , (22) 

where /o is the initial population, which can of course 
be taken to be zero. However, notice that the popu- 
lation reaches the saturation limit fs independently 
of the initial value /q. This contrasts with the solu- 
tion of the GL model 7 , which requires a non-zero 
value /o or otherwise the solution is the trivial solu- 
tion f{t) = 0. 

4.2 Dynamic environment 

We consider the case in which the environment can 
change due to external influences. For us this means 
that the parameters 7 and fs that appear in the 
model equation, both of which depend on the state of 
the environment, change with time. In the absence 
of a dynamical theory of the interactions of the cells. 
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all we can do is promote 7 and fs to be functions of 4.3 Examples 
time. Therefore, our basic equation becomes 



dt 



-lit) [I-fs{t)] 



(23) 



This equation is conveniently solved by the Green 
function method, 

/•OO 

/(<) = / dt'G{t,t')-i{t')f,{t') + hit) , (24) 
Jo 



where G satisfies 



rir' 



(25) 



fh is a solution to the homogeneous equation such 
that / satisfies the initial condition /(O) = /q. In 
Eq. (|25|l . 6{x) stands for the Dirac delta function. A 
suitable Green function for Eq. (|23|l is 

G{t,t')^e{t-t')g{t,t'), (26) 

where g{t, t') is the solution to the homogeneous 
equation satisfying the condition g{t',t') = 1, and 
6{x) is the unit step function. The function g{t, t') is 
then uniquely determined as 

g(i,0^e-/>*"''(*"\ (27) 
and therefore the solution for f{t) is given by 



fit) = / dt' git, t')-tit')fS') + /o git. 0) . (28) 
Jo 

Needless to say, if 7 and fg are assumed to be con- 
stant, then Eq. H28|l reduces to the solution given in 
Eq. H22() . But for any imaginable functions 7(t) and 
fs it) that could be used to parametrize the changing 
environment, Eq. H28|l readily provides the complete 
solution in terms of two integrals. 

A particularly simple form of the solution is ob- 
tained in the case in which fs is a constant. In this 
case, noting from Eq. (|27() that 



'^^|^=7(t')5(M'), 



dt' 



then Eq. (gHJl yields 

/(t) = /,[l-5(t,0)] + /o5(t,0), (30) 
where we have used git, t) — \. 



In order to illustrate some general features of the so- 
lution, we will consider below various specific cases. 

4.3.1 Example 1 

Let us assume that fs is a constant, while 7 varies as 
some (integer) power of t; i.e.. 



fs — constant 

7 = af\ 



(31) 



where a is positive constant and n is a positive inte- 
ger. First, from Eq. (|77|l . 



5(M') = e--(*"^^-'"^0, 
where we have defined 



70 = 



71+ I 



(32) 



(33) 



for simplicity of the notation. The solution obtained 
from Eq. (|^ is then 



fit) = fs 



1 - e 



+ foe 



'lot" 



(34) 



In Fig. n we plot the function f/fs, in arbitrary 
time units (i.e., setting 70 — 1), for the values of the 
exponent n = 1, 2, 3, 4, 5, and taking /o — 0. 

4.3.2 Example 2 

As before, we assume that fs is a constant, but now 
take 7 as a combination of two monomials; i.e.. 



fs = constant 
7 = at" + bt"" . 



(35) 



where n, m are positive integers and a, b are positive 
constants. Following the same steps as above, the 
obtained from Eq. H3()(l is given by 



(29) fit)^f 



,-(7ot"+'+7ot'"+') 



(36) 

where 70 is defined as in Eq. (|33|l and in analogous 
fashion 

h 

(37) 



7^^ 



TO + 1 



5 




Figure 1: Plot of f{t)/fs with f{t) given by Eq. 1^ . 
with the initial condition /o = 0, for 7 — at" with 
n — \,2, 3, 4, 5, in arbitrary time units. 



Figure 2: Plot of f{t)/fs with /(t) given by Eq. 
with the initial condition /o = 0, in arbitrary time 
units, for various values of n, m and the ratio r = 

7o/7o- 
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In Fig. 12 we plot the function ///s, for various values 
of n, TO and the ratio r = 70/70 • 

We stress that, apart from the initial and satura- 
tion values /o and fs, respectively, the only unknown 
and ajustable parameters of these models are the con- 
stants that appear in the ansatz for the growth rate 
function 7(t), e.g. the constants a and 6 defined in 
Eqs. (EB and 123). 

4.4 Generalizations 

By inspection, Eq. (|36|l can be generalized in an obvi- 
ous way to the case in which 7(f) is a power series in t. 
Although the models considered above, with fs taken 
to be a constant, already provide a flexible and rich 
structure as far as their ability to fit the phenomen- 
logical data is concerned, it is possible to consider the 
cases in which fs is a function of i. In these cases it is 
not possible to give a closed expression for the solu- 
tion given in Eq. (|28|l . in general. However, in other 
analogous physical problems where similar situations 
arise|S], very effective approximation methods have 
been used which could be employed in these cases as 
well. 

5 Conclusions and Outlook 

In this article we have presented a model for tumor 
growth, which is based on physical principles on one 
hand, together with plausible physical assumptions 
and idealizations on the other. As shown in Section 
01 in the simplest version of the model, in which the 
pairwise interactions between the cells are neglected, 
the growth equation is linear and a complete solution 
can be readily given. Moreover, the solutions were 
explicitly given for a few sample cases, which exhibit 
the known characteristic features of tumor growth. 
Thus the approach that we have followed is fruitful 
in several ways. Firstly, the model contains a min- 
imal number of parameters, which have a concrete 
and well-defined meaning, and arc in principle deter- 
mined and calculable if the interactions that govern 
the underlying microscopic mechanism of growth arc 
known. Secondly, by the same token, the use of this 
model should in turn shed light on such mechanisms. 



thereby providing a firm basis for pursuing this line 
of work. Thirdly, the model can be extended beyond 
the linear approximation that we have used if the 
pairwise interactions are believed to be important in 
a particular system and the quadratic terms in Eq. 
(|^ should be taken into account. Our approach pro- 
vides a framework for taking into account such higher 
order terms in a systematic fashion. Lastly, while in 
this paper we have restricted ourselves to treat the 
growth of an untreated tumor, our work paves the 
way for applying similar principles and ideas to in- 
clude the efi'ects of therapy. Work along these lines 
is in progress. 
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